EXPLICIT INVERSIONS OF CERTAIN MATRICES I 



RUIMING ZHANG 



Abstract. In this note, we demonstrate a method to invert some Hankel 
matrices explicitly by using the kernel polynomials for the related classical 
orthogonal polynomials. 



1. Introduction 

In the theory of orthogonal polynomials, We could calculate the determinants 
of some Hankel matrices once we know the three term recurrence relation for the 
associated orthogonal polynomials and vice versa. It is well-known that the kernel 
polynomials of the orthogonal polynomials encodes important information about 
the associated Hankel matrices. These matrices are generalizations of the Hilbert 
matrices. In this note we present a method to invert some Hankel matrices associ- 
ated with classical polynomials by using the kernel polynomials. 

The following theorem very is a well known fact from the theory of orthogonal 
polynomials: 



Theorem 1.1. Given a probability measure fx on K with a support of infinite many 
points. Let us consider the Hilbert space of [i-measurable functions 



(1.1) 



X 



f(x)\ / \f(x)\ 2 d^x) < oo 



with the inner product defined as 

(1.2) (/,«?) := [ f(x)^)d^(x), f,geX. 



Assume that {w n (x)}'^' =0 is a sequence of linearly independent functions in X with 
wq(x) = 1. Let 



(1.3) 



(1.4) 



ano 



Wj(x)wk(x)d[i(x)i j,k = 0,l. 



n„ 



aio an 
\ a n0 ot n i 



Qnn J 



ne NU{0}. 
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(1.5) 



A„:=detLT„, n <E N U {0} . 



Then the n-th orthonormal function with positive coefficient in w n {x) is given by 
the formula 



( «00 "01 



(1.6) Pn(x) 



: (let 



&02 



a w an «12 



a„ a„i a n2 
\ w (x) wi(x) w 2 {x) 

for n G N , with 

(1.7) p (x) = w (x) = 1. 

Furthermore, the coefficient of p n (x) in w n (x) is 



Otln 



w n (x) J 



(1.8) 



7n : = 



A„_i 



Proof. The proof is the same as for the case w n (x) = x n , which could be found in 
any orthogonal polynomials textbooks such as [2]. □ 

Corollary 1.2. For n € N ; we have 

n 



(1.9) 



k=l 



In 



Proof. This is a trivial consequence of (|1.7p and (|1.8j) . 
Lemma 1.3. Let 



□ 



(1.10) 



y) :=J2Mx)Pk{y), neNU{0}. 



Then, for any it(x) in the linear span of {wkix)}^ , we have 



(1.11) / Tr(x)kn(x,y)dfi(x) = n(y). 

Proof. To see (|l.lip . just expand 7r(a:) in Pk{x), k = 0, 1 . . . , n. □ 
Lemma 1.4. For each n S NU{0} ; the function k n {x, y) satisfying {LTT| is unique. 
Proof. Suppose there are two such functions h n (x,y) and k n {x,y) , then, 

(1.12) < ||^^(-,2/) — fc^(-,2/)|| 2 

={h n {-, y) ~ k n (-,y),h n (-,y) - k„(-,y)) 

=(h n {-,y) - k n (-,y),h n (-,y)) - (h n (-,y) - k n (-,y), k n (-,y)) 

=0, 

which is a contradiction. □ 
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Lemma 1.5. Let 

(1.13) (fok)o<j,k<n ■ = n-\ n e N U {0} . 

Then, 

n 

(1.14) k n (x,y) = Pj k Wj(y)w k (x). 
Proof. Let 

n 

(1.15) f(x) = }^UkWk(x), 

k=0 

then, 

n 

(i-i6) (/(•), £ PjkwffiM-)) 

j,k=0 

n 

= ^2 U m(Wm(-),k(-,y)) 
m=0 
n n 

= E Mm E PjkWj(y){w m ,w k ) 

rn=0 j,fc=0 
n n n 

m=0 j=0 k=0 



By Lemma [Ql we have 

(1.17) k n (x,y)= Pj k Wj(y)wk{x). 

j,k=0 

Corollary 1.6. TTie kernel in Lemma f"OI is also given by 



(1.18) K{x,y) = -— det 



/ 1 

1 aoo "oi 

\ w n (x) a„ a„i 



/or 7i E NU{0}. 
Proof. Since 

n 

(1.19) k n (x,y) = ^2 P]kWj{y)wk{x) 

j,k=0 

with 
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Then, 
(1.21) 



n„(fc,j) n„(fc,j) 



Jk detn„ A„ 
where n n (fc,j) is the (k,j)-th co-factor. Therefore, 



(1.22) 

It is clear that 
(1.23) 



I « 

k n {x,y) = — ^2 u n(k,j)wj{y)w k (x). 



j,k=0 



^2 Tl n(k,j)w j (y)w k (x) 
j,k=o 



(W(y)) 

w(x) n„ 



by direct determination expansion, which is 



(1.24) 
where 

(1.25) 

and 
(1.26) 



1 



W(x) 



(W(y)) 

w(x) n„ 



( 1 \ 

\ w n (x) J 



(W(y)) = ( 1, Wl (y) , w n (y) ) 



□ 



Lemma \l .51 enables us to compute the inverse the Gram matrix in terms of the 
orthonormal functions {pn(x)} n — - 

Corollary 1.7. Assume that {w„(x)}^ , {p n (a;)}^ and II„ = {otjk)o<j,k<n as 
in Theorem \l.l[ Suppose we have two families of linear junctionals {wfc}^L an d 
{ v k}kLo over li near space generated by {wn( x )}™=o w ith 



(1.27) 
and 



Uj(w k ) = 5jk, 



(1.28) 

for j, k = 0, 1, . 



Vj(w k ) = 8jk 



Then, 



(1.29) 

where 
(1.30) 



/%* = U k(Pm{x))Vj(p m (y)), 



'jk)0<j,k<n 



n: 
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Proof. From Lemma 11.51 we have 

n n 

(1-31) ^2 f3 jk Wj(y)w k (x) = ^2 Pm(y)Pm(x). 

j,k=0 m=0 

Then we apply the functional Uj and v k both sides of the above equation, the claim 
of the corollary follows. □ 

2. Main Results 
2.1. Preliminaries. The Euler's T(z) is defined as [I] 



(2.1) 



1 



T(z) 



n 1 



3=1 



1 



z e 



For a, oi, a r £ C, the shifted factorials are defined as 
T{a + n) 



(2.2) (a), 



r(o) 



(ai,...,a r )„ := JJ(oj)n) neZ,reN. 

3 = 1 



The generalized hypergeometric series r F s with parameters {ai, a r } and {&i, 6 S } 
is formally defined by 



(2.3) 



at, a,2, a r 



E 



[at, ■ a r ) n z n 



^ (oi, ...,a s ) n n\ 



bt,b 2 , —,b s 
The Barnes G-function is defined as 

oo 

(2.4) G(z) : = (2^/ 2 e -[^+D+^ 2 ]/2 TT (1 + £)» e -.+*V(2») 

- L - L n 

n=l 

where 
(2.5) 



7 



/ " 1 \ 
= lim I > - — In n . 

n— >oo \ ^ — ' fc / 

\fe=l / 



The Barnes G-function is an entire function with the property 
(2.6) G(z + 1) =T(z)G(z), 



(2.7) 

and 

(2.8) 



fe=0 



Gin) 



n = 0,-l,-2,. 



In following, in all the cases except the last one we use functional 



(2.9) 

and for the last case we use 



Ui(p(x)) = Vi(p(x)) = i 



d l p(x) 



(2.10) 



i(p{x)) = Vi( P (x)) = - 



dx l 



d l p(x) 



x=0 



dx l 



where p(x) is a polynomial in variable x. 
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2.2. The Hermite Polynomials. The Hermite polynomials {H n (x)}^°_ are 
fined as pQ 



(2.11) 

for n > and 
(2.12) 

They satisfy 
(2.13) 



n n I 1 

2 ' 2 ' 2 



ff„(a;) = (2a;)" 2 F 

= 0. 

DH n {x)=2nH n _ 1 (x), neNU{0}. 



Hermite polynomials satisfies 



(2.14) / H n (x)H m (x)exp{~x 2 )dx = 2 n n!M 

7k 



for n, m = 0, 1, .... 

Thus, the orthonormal polynomials 

(2.15) h n {x) : = 
have leading coefficients 

(2.16) 7n = 
Clearly 

/oo 
y n e- y2 dy 
-oo 

and 



_Hn(x)_ 



l + (-l) n T fn+l 



(2.18) aij 
for i,j = 0, 1, ...,n. Thus, 
(2.19) 
or 



l + (-l)' +J r + j + 1 



detl 1 + ( : 1)<+ W' + J " + 1 



j,fe=0 



2 5— 7T— ]__[ fc! 

fc=0 



deti 1+( : 1) ' + W* +J+1 



(2.20) 



for n = 0, 1, .... 

The (i,i)-th entry of If- 1 = (/3 jfc )" fc =o is 



= 2" Ziil ^ til 7r 2! * i G'(n + 2) 



j,fc=0 



1 n 1 



fc!2 /c v / 7T 

fc-max(ij) 



£C = 



v=o 
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or 



(2.22) ftj = ^ ± ^ {(')&-.(»)} {0)^(0)} 
for i. j = 0, 1, n. 

Theorem 2.1. For n € N U {0}, i/ie matrix 

(2.23) ^ — r' J 



2^ 



0<i,j<n 



has the determinant 

(2.24) de «(i±i±Cir(i±i±^y =2-^ G (n + 2) , 



ij=0 



and iis inverse matrix is 



(2-25) ( £ L ^ 

, /c— max(i,j) 



0<i ,j ' <n 



2.2.1. T/ie Laguerre Polynomials. The Laguerre polynomials {£n( x )}^Lo ma y ^e 
defined as [I] 

(2.26) L-(«) = ( ~^ ;x 

for n > 0, and we assume that 

(2.27) £-i(z) = 0. 
We also have 

(2.28) d J^ = - K +\ [x) 
and 

(2.29) I" L a m {x)L a n (x)x a e- x dx = ^ + " + ^ S mn , 



o 



for a > — 1 and n, m = 0, 1, .... The orthonormal polynomials 



(2.30) ^ (8) = ( - ir ^ _^_ £ g (!B) , 

have the leading coefficients 

(2-31) 7n = , , , 1 = 

% /n!r(a + n + 1) 

for n = 0, 1, .... Clearly 

i>00 

(2.32) / x n+a e- x dx = T(a + n + l), 

Jo 

and 
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(2.33) 

for i,j = 0, 1, ...,n, Then, 



Oij — r(a + i + j + 1) 



det (T(a + i + J + l))" fe=0 = I] { fc!r ( a + fc + !)> ' 

fc=0 

G(n + 2)G(a + n + 2) 



(2.34) 
or 

(2.35) det(r(a + i + j + l))£ fc=0 

Let n- 1 = (Pjk)%k= > then ' 



G(a + 1) 



n 



it! 



— r(a + fc+l) 



efL£(x) 



i=0 



" i.! r 

E r^TlTT) K^^U [(-i)^B(v) 

nax(jj') 



or 



(2.37) 



' iij! ^ r(a + fc + 1) 

/c-max(tj') 



for j, fe = 0, 1, n. 

Theorem 2.2. For n = 0, 1, the matrix 
(2-38) ((a+l) i+j ) <^<„ 
/ias determinant 

G(n + 2)G(a + n + 2) 



(2.39) 
and inverse 

(2.40) 



det((a + l)i +J )". = 



G(a + l)r(a + 



En (a+l) fc ffc^ /fc 

fc=max(i,j) fc! 



(-l)<+*(a +!),•(* +1), i 



0<i,j<n 

2.2.2. The Ultraspherical Polynomials. The Ultraspherical polynomials (or Genen- 
bauer polynomials) {C^(x)} n _ are defined as, [1] 



(2.41) 



(2A) n 



^1 



-n, 2A + n 1 - a; 



for n > 0, and we assume that 

(2.42) C\{x)=0. 

We also have 



(2.43) 



(2.44) 



2\\- 



?rr(2A + n) 



2 2A - 1 n!(A + n)[r(A)] 2 



for A > — i and n, m = 0, 1, .... The orthonormal polynomials 
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(2.45) Pn (x) 
have leading coefficients 



(2.46) 

It is clear that 



7n 



' 2 2A ~ 1 n!(A + n)[r(A)] 2 A 

7rr(2A + n) ° nW 



'(A + n)2 2A + 2 "- 1 r(A + n) 2 



7m!r(2A + n) 



(2.47) 
and 



(2.48) 



for i, j = 0, 1, ...,n, where B(p,q) is the beta integral 



(2.49) 



B(p,q)= [ xP-^l-x^dx, R(p),R(g) > 0. 



Then, 
(2.50) 



(l ,. 1 ii±i^ s ri±4±i,A + ii)=n 



7rfc!r(2A + /j) 



fc=0 



(A + fc)2 2A + 2fe -!r(A + fc) 2 



7r" +1 G(n + 2) G(2A + n+ l)G(A) 2 



or 

(2.51) 

i , i i +(- i ) i+i o/ i +j+ i x , r 

"■" ' ' 2 V 2 ' 2JJ 2(«+ 1 )(«+ 2A -i)(A) rl+1 G(2A)G(A + n + l) 2 - 

The (i, j)-th entry of the inverse matrix II" 1 = (/3jfc)™fc = o ls 
(2.52) 



1_ ^ 2 2A - 1 fc!(A + fc)[r(A)] 2 [d 4 G A (x) 



= 7T7T E 



fc-max(ij') 

2 2A -![r(A)] 2 



7iT(2A + fc) 



fix 1 



x=0 



2/=0 



i\j\Tr ^ . r(2A + fc) 

/c-max(i.j) 



j/=0 



2 t+J (A)i(A)jT(A) ^ fc!(A + fc)G^(0)C^(0) 



(2A) fc 



for i,j = 0, 1, ...,n. 
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Theorem 2.3. Let (aij)o<i,j<n be the matrix with entries 

l + (-l) t+j + j + 1 , 1 
(2-54) a ij = '-^-^Bl^—L ,A + - 

for i, j = 0, 1, n, then, 
(2.55) 

+ / f + j + 1 m " 7rfc!r(2A + fc) 

d6t v 2 * \~^2~' X + 2 J J " 11 (A + fc)^+2 fe -i r(A + fc)2 • 

07' 

(2.56) 

;); , [ " l + (-l) i+J g + 1 A _ ^ l+1 G(n + 2) G(2A + n+l)G(A) 2 



2 V 2 2JJ 2("+ 1 )(™+ 2A - 1 )(A)„+i G(2A)G(A + n + l) 2 

T/ie inverse matrix (/3ij)o<i,j<n has entries 

/ori,j = 0, l,...,n. 

2.2.3. The Jacobi Polynomials. The Jacobi polynomials { P^'^lx) \ may be 

I J 71=0 

denned as [HE] 

(2.58) P^(x) = ^±^ 2 F 1 ( - n; ^J+ /3 + 1 ;^) 
for n > 0, and 

(2.59) P ( _ a { P) (x)=0. 
We also have 

(2.60) dP - a 2 {x) = n+ ^ +1 er"H, 



and 



(2.61) | ll" ^P^^M^dx = ft n 5 mn 

for a,/3 > —1 and n, m = 0, 1, ... with 

(2.62) w(x) : = (l~x) a (l + xf, 
and 

2 a+/3+lp( a + „ + + n + 1) 



(2 ' 63) hn '' {2n + a + f3+l)T(a + f3 + n+l)n\' 

The orthonormal polynomials 



(2.64) p n {x) 



l (2n + a + [3 + l)T{a + (3 + n + jjnj (a|g) , , 
V 2«+' 3 + 1 r(a + n + l)r(/3 + n + 1) " W 
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with leading coefficients 

(Z.DO) 'Tn - . 

^/nlTrl^a + /? + n + l)r(a + n + l)T([3 + n + T) 
for n = 0, 1, .... The moments of the Jacobi measure are 



(2 - 66) M " = (-i)»r(a + /?+l) 2 ^ 



-n,(3+l 
a + (3+ 1 



for n = 0, 1, .... Then (i, j)-th entry of the Hankel matrix H n = (ctjk)™ k=0 is 

-i-j,0+l 



(2.67) 



(-l)'+JT( a + l)rQ3+l) 

ij - ^ w . „ . ,x 2fl 



2-«-/?-ir(a + /3 + 1) 
for i, j = 0, 1, n. Then, 



;2 



(2.68) 



detn„ - 



G(n + 2)G (^+|+i) 2 G (^±^V 



2 y 



(^±|±i) n+iG (^±|±3 + n) 2 G(^ + ny 
G(a + /3 + n + 2)G(a + n + 2)G((3 + n + 2) 



G(a + /? + l)G(a + l)G(/? + l) 
and the (i, j)-th entry of the inverse matrix II" 1 is 
(2.69) 

a J_ A (2fc + a + /? + l)r(a + g + fc + l)fc! 

2 a +f>+ 1 T{a + k + 1)T(J3 + k + 1) 









dx l 


x=0 


dy j 



1 ™ 
"7T7T X] 



(2fc + a + (3 + l)r(a + /3 + k + l)k\ 



^ 2 a +P+ 1 T(a + k + l)T([3 + k+ 1) 

fc-max(i.j) 



(fc + a + P + 1)^ 



(a+i,/3+i) 



(X) 



x=0 



2i 



{a+3,0+3) ^ 



7=0 



or 



(2.70) fa = J2 

fc=max(i,j) 



(2fc + a + /? + l)r(q + P + fc + l)fc! 
r(a + + 1)T{(3 + k + 1) 



\ {k + a + (i+ l) i P fc ( y /?+z) (0)}{(fc + a + j3 + l)j^y^(0) 

X | j!j! 2 a+/3+i+j+l 

for i, j = 0, 1, ...,n. 

Theorem 2.4. For n = 0, 1, ifte matrix 

-i- 3,0+1 . 



(2.71) 



•A 



a + /3+ 1 



0<i,j<n 
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has the determinant 



(2.72) det( 2 F l( a+ y +1 ;2 



j,k=0 
n+l 



i- 3,0 + 1 
oi + + 1 

T(a + + l) 2 ~(2"+2/3+n+l) 7r 

r(a + / 9+ l)r(a + /9+ 1) 
G („ + 2 )G (^±f±i) 2 G f 2 



G (£±|±3 +71 ) 2 G (£±|±4 +n N2 

G(a + + n + 2)G(a + n + 2)G(0 + n + 2) 
/a+£+i\ G ( a + /3+ i) G ( a + i)G(/3+l) : 

V z /n+l 



n+l 

and iis inverse matrix (0ij) o<i i<n with 



(2.73) /3y = J2 



0<i,j<n 

k\(2k + a + 0+l)(a + + l) k 



(-2)*+ii!j!(a + l) fc (/3 + l) fc 

K-max(ij) 

xj (£; + « + /? + l) 4 Fi_t 4,/3+l) (0)}{(fc + a + + l)A-f f5+1 \Q)} ■ 
fori,] =0,l,...,n. 

If we take the polynomial sequence 
(2.74) w n (x) = (x-l) n 

for n = 0, 1, and the linear functional defined in <|2. 10|) . Then the (i, j)-th entry 
of n n = (a jk )l k=0 is 



el 

(2.75) a# = I (x- iy +1 w{x)dx 



_ 2°+^+ i +J+ 1 r(a + i + j + l)r(/3+l) 
( ?bj (_i)i+ir(a + /9 + i + i + 2) 

for i,j = 0, 1, ...,n, and its determinant is given by 

G(a +n+ 1)G(0 + n + l)G{a + + n + 1) 



(2.77) detn„ 



2 (n+a+0)(n+l)G( a + 1)G(0 + 1)G(q + + 1) 

\ 2 
a±/3+2\ 



7r "+ 1 G(n + 2)G(^^±§±±) g( 



(s±|±i) G(s^±i+»+l) G(^±f±2+n+l)' 
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and the it inverse matrix has entries 
(2.78) 

(2k + a + /3 + l)T(a + + k + l)k\ 



i 



<-^> 2 a +P+ 1 T(a + k + l)T(/3 + k + 1) 

fc=max(i,j) 









dx l 


x = l 


dy j 



E 



fc=max(i,j) 



(2fc + a + (3 + l)r(a + /? + fc + l)fc! 

2 a+/3+i r ( a + fc + i)r(/? + jfc + 1) 



{k + a + + l) t PtV 4)+l \x) 
i\2* 



J x=l L 



./'!-" 



Jy=l 



or 



^ (2fc + a + /3+l)r( a + /3 + fc + l)fc! 
( j / * y ~ ^ , r(a + * + l)r(/J + * + l) 

fc-max(t,j) 

r (fc + a + a + i),p^ +i) (i)}{(fc + « + + ' 



for i,j = 0, 1, n. Since 



(2.80) 
then 



p(a,/3) m _ ( a + 

n lJ ~ rc! ' 



(2.81) 



Tja + /? + l)(g + /? + l^g + /? + l)j 

2 a+/3+i+j+l( Q! + ^.(q, + !) jr ( a + + 
(2fc + a + /3+l)(a + l) fe ' 



fc!(a + /?+l) fc (/3 + l)fc 



x E 

fc-max(ij) 

(*](«. f f / - ! )/.(n - / + ., + i ),, 



for i,j = 0, 1, ...,n. Therefore, we have proved the following: 
Theorem 2.5. For n = 0, 1, the determinant of the matrix 

(2.82) ( {a + 1) ^ 



( a + P + 2 )i+iJo<ij< n 
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(2.83) det 



(a + /3 + 2) i+j J ij=0 



G(^±|±i) g(* 



2 



G(a + l)G(/3 + l)G(a + /3+l) 
7rr(a + /3 + 2) 



n+l 



2 2r l +2a+2/3+i r ( Q , + + !) 

G{n + 2)G{a + n+ l)G(/3 + n + l)G(a + /3 + n + 1) 



{ 2± § ±1 ) n+1 G( 2± § ±1 +n+l) G(?±§^+n+l) 
and its inverse matrix (jij) 0<i j <n has elements 



2 " 



(2.84) 7y 



(a + l)i(a + l)j(a + /3+l) 



fc=max(i,j) 



(2fc + a + /3+l)(g + l) fe 
*% + /?+ l)k(/?+l)k 



k\ (k 

' ' . )(a + p + i + l) k (a + (3 + j + l) k 



fori, j = 0,l,...,n. 
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